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The connection between work and changes in the Hamiltonian for a system with a time-dependent 
Hamiltonian has recently been called into question, casting doubt on the usefulness of the Jarzynski 
equality for calculating free energy changes. In this paper, we discuss the relationship between two 
possible definitions of free energy and show how some recent disagreements regarding the applica- 
bility of the Jarzynski equality are the result of different authors using different definitions of free 
energy. Finally, in light of the recently raised doubts, we explicitly demonstrate that it is indeed 
possible to obtain physically preprintrelevant free energy profiles from molecular pulling experiments 
by using the Jarzynski equality and the results of Hummer and Szabo. 
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Single molecule experiments, such as the stretching of 
a polymer molecule using an atomic force microscope or 
laser tweezers, have become common in the last decade 
P, The goal is often the determination of the free 
energy surface along some coordinate of the molecular 
potential energy surface. In order to proceed, one in- 
vokes the Jarzynski equality using the extension derived 
by Hummer and Szabo [3]. 

Although there has been some controversy about these 
theoretical advances, it is fair to say that their use in 
interpreting nano-scale single molecule experiments is 
widespread. Thus any question that they may be fun- 
damentally in error must be carefully examined. 

Recently, questions have been raised about the con- 
nection between work and changes in the Hamiltonian 
for a system with a time-dependent Hamiltonian, cast- 
ing doubt on the applicability of the Jarzynski equality 
for computing free energy changes [1]. Here, we discuss 
these questions and show that the Jarzynski equality can 
be usefully applied to determine physically relevant free 
energy changes. 

Consider a system with Hamiltonian 7io(a;), where x 
represents the microstate of the system, and suppose that 
this system is subject to a time-dependent force f{t) act- 
ing along some coordinate z{x). From the perspective 
of classical mechanics, we have two options for treating 
the force. We may consider it as an external force not 
included in the Hamiltonian of the system and study the 
evolution of a system governed by TLoix) under the effect 
of the external force f{t) acting along z{x). Alterna- 
tively, we may include the force in the Hamiltonian of 
the system and study the evolution of a system governed 
hy n{x,t) =no{x) - z{x)f{t). 

In the first case, we are considering a time-independent 
Hamiltonian under the effect of an external force f{t). 
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According to classical mechanics, the work done by the 
external force up to time r is 



W{t) = / f[t)d{z[x{t)]} 

Jt=Q 



(1) 



and we have the usual result that the work done on the 
system equals its energy change. 



no[x{t2)] - Ho[x{ti)] = W{t2) - W{ti). 



(2) 



The free energy change appropriate for this first descrip- 
tion of the system is 
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(3) 



{ksT = 1 throughout). 

In the second case, we consider the time-dependent 
Hamiltonian TL{x,t) — Ti.o{x) — z{x)f{t). In this descrip- 
tion of the system, f(t) is an internal force and there 
should be no expectation that the work done by /(t) 
equals the change in energy of the system. Here we con- 
sider the thermodynamic work, 



Wtir) 



dt 



m 



(4) 



which by definition equals the change in energy of the 
system. The appropriate free energy change to consider 
for this description of the system is 



AGt(T) = Gt(r)-Gt(0) = -log 



(5) 



We note, as suggested by Vilar and Rubi (VR) [J], that 
this second description of the system is not unique - 
adding a term g{t) to the Hamiltonian has no effect on 
the dynamics of the system but changes the values of Wt 
and AGi. 

In considering the effect of a force / on a harmonic 
spring of force constant k, VR describe the system ac- 
cording to the first picture and obtain AG ^ W ~ /^/2fc 
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while Horowitz and Jarzynski (HJ) use the second pic- 
ture and obtain AG* = Wt = -f/2k 0, Q. Both of 
these results are correct in their respective descriptions, 
and mean different things. In particular, VR are describ- 
ing the free energy change associated with changing the 
length of the spring in the absence of an external force; 
the force is only a tool used to measure the free energy 
profile of the free spring. Meanwhile, HJ are describ- 
ing the free energy change of the combined force-spring 
system as a function of the force. 

The Jarzynski equality is framed in the second of our 
descriptions and expresses a relation between Wt and 
AGt 



-AGt _ /^-Wt 



The validity of this expression is not in question - only its 
utility in describing free energy changes in a system. VR 
point out that AG* depends on the arbitrary choice of 
g(t) in the Hamiltonian and leads to arbitrary free energy 
changes. If all that can be extracted from the Jarzynski 
equality is this arbitrary AGf, then the Jarzynski equal- 
ity seems to be of little use. We shall show, however, that 
this is not the case. 

Consider a single- molecule pulling experiment, for 
which the Jarzynski equality has frequently been applied 
P, 0| • In studying the unfolding of a biomolecule, one is 
often interested in the free energy profile G{z) as a func- 
tion of end-to-end distance z. We could map the free 
energy by reversibly pulling the ends of the molecule and 
measuring the work exerted by the external force as a 
function of z. This is of course the classic method and 
corresponds to VR's analysis of the harmonic spring. 

We could also try to get the free energy profile using 
the Jarzynski equality. Direct application of the Jarzyn- 
ski equality to yield AGi gives the free energy difference 
between the free molecule and the molecule with a cer- 
tain force applied to it. This is not in itself a particularly 
useful quantity and is not the free energy profile. Hum- 
mer and Szabo have, however, shown how to obtain free 
energy profiles from single-molecule pulling experiments 

Consider an unperturbed system described by a Hamil- 
tonian Ti.(){x). When a time-dependent perturbation is 
applied along some coordinate z{x), we write the new 
Hamiltonian as Ti.o{x) + Ti.'{z,t). Hummer and Szabo 
have shown that the unperturbed free energy profile 
along coordinate z can then be reconstructed as 

G{zo) = - \og{S[zit) - zo]e'^'+^'(^'*^ (7) 

where the average is over all trajectories of the system in 
the presence of the perturbation Q . 

We now apply this result to a macroscopic, determin- 
istic spring and show how the Jarzynski equality can be 
used to calculate G{z), thereby reconcihng the results of 
HJ and those of VR. In this case, there is only one de- 
gree of freedom so the microstate x is simply the length 



of the spring. Our model is Hoix) = + kx^ /2 and 
TL{x, t) = p"^ + kx'^/2 — f{t)x, where f{t) switches from 
to /o over time < i < r. The pulling process is finished 
at t = T, at which time /(t) — /q and from classical 
mechanics x(t) = fo/k. 

Inserting this into Hummer and Szabo's result, 

G{z) = - \og{S[x{T) - z]e-^''f'>^^^^), (8) 

where we have used the fact that in this case our coordi- 
nate of interest z{x) is just x. 

In this deterministic case only z — fo/k contributes to 
the average and we have 



(6) G{z = fo/k) = - log(e-'^'-^o/'=) = - log \{e-^')e-^o/k 

(9) 



For reversible pulling, HJ calculate that the work distri- 
bution is sharply peaked at Wt = — /Q/2fc, so we finally 
obtain 



G(z = /o/fc) = -log 



./o/2fcp-/o7fc 



= +/oV2fc, (10) 



which agrees with VR and is the expected result for the 
free energy profile of a Hookean spring. 

Consider the effect of adding an arbitrary g{t) to the 
general Hamiltonian. The effect on AGt is easily seen 
from Eq. (O to be 



AG— (r) = AGt(r) + [5(r)-g(0)]. 



(11) 



Since the term g(t) redefines the zero of energy at each 
point in time, it is expected that AGt(T) will be affected 
as it is comparing free energies at two different times. 
Before ascribing a physical interpretation to AGt , it must 
be corrected by subtracting this arbitrary change in the 
zero of energy. 

We now examine the effect of an arbitrary g{t) on the 
free energy profile G(z) computed via Jarzynski's equal- 
ity. We then have 



'ix,t) = n'{x,t)+git) 



y^new 



Wt+g{r)~g{0) 



(12) 
(13) 



and Eq. ([7|) becomes 

G"^'"(zo) = - log{S[z{t) - zo]e-'^'+^'(^'*'+»(°^), (14) 



which can be simplified to 



G'-™(zo) = G(zo)-5(0). 



(15) 



So adding a time-dependent term g{t) shifts the overall 
free energy profile G(z) by an additive constant, but has 
no effect on relative free energies. 

We can consider the same situation from the perspec- 
tive of thermodynamics in one dimension. The internal 
energy of the system is given hy U = TLq and its enthalpy 
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by the Legendre transform Ho — fz. We can then define 
two free energies, G ^ U - ST and Gt = U - fz - ST. 

The method of VR is constructed to measure G as a 
function of position, G(z), while HJ arc calculating Gt as 
a function of /, Gt{f ). As long as the fluctuations in x 
are small at a given / (ie, we are in the thermodynamic 
limit), we can use the simple relation G{z) = Gt{f) + fz 
to convert between the two quantities. Outside of this 
limit, there is not a simple relation between the two quan- 
tities but the method of Hummer and Szabo discussed 
above can be used to reconstruct G{z) from pulling ex- 
periments. 

In conclusion, we have shown that by properly apply- 
ing the Jarzynski equality, the textbook result for the free 
energy profile of a spring is correctly recovered. More im- 
portantly in light of recent doubts, we have reaffirmed the 
applicability of the Jarzynski equality to the analysis of 
single-molecule pulling data. 
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